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ABSTRACT 

Aims. To derive a new map of the interstellar extinction near the Galactic Center (GC) extending to much higher values of Ay than 
previously available. To use the results obtained to better characterise the long-period variable star population of the region. 
Methods. We take the Spitzer IRAC catalogue of GC point sources (Ramirez et al. 2008) and combine it with new isochrones (Marigo 
et al. 2008) to derive extinctions based on photometry of red giants and asymptotic giant branch (AGB) stars. We apply it to deredden 
the LPVs found by Glass et al. (2001) (Glass-LPVs) near the GC. We make period-magnitude diagrams and compare them to those 
from other regions of different metallicity. 

Results. Our new extinction map of the GC region covers 2.0 ° x 1.4 ° (280 x 200 pc at a distance of 8 kpc). The Glass-LPVs follow 
well-defined period-luminosity relations (PL) in the IRAC filter bands at 3.6, 4.5, 5.8, and 8.0 /jm. The period-luminosity relations are 
similar to those in the Large Magellanic Cloud, suggesting that the PL relation in the IRAC bands is universal. We use ISOGAL data 
to derive mass-loss rates and find for the Glass-LPV sample some correlation between mass-loss and pulsation period, as expected 
theoretically. Theoretical isochrones for a grid of different metallicities and ages are able to reproduce this relation. The GC has 
an excess of high luminosity and long period LPVs compared to the Bulge, which supports previous suggestions that it contains a 
younger stellar population. 

Key words. ISM: dust, extinction - Galaxy: stellar content - Infrared: stars 



1. Introduction 

The variable stars of the Galactic Center (GC) region are of 
great interest, both as population and distance indicators. Long- 
period, large-amplitude variables situated on the asymptotic gi- 
ant branch (AGB), comprising Miras and OH/IR stars, are the 
easiest objects to detect thanks to their high luminosities. Thus, 
they are one of the few stellar populations that can be observed 
in their entirety towards the GC. As is well-known, in the inner- 
most parts of the Galaxy, surveys are hampered by high inter- 
stellar extinction (e.g. Frogel et al. 119991 Schultheis et al. 1999); 
observations must therefore be carried out in the infrared where 

can be as low as 0.04 Ay for A/um < A < 8jum (Indebetouw 
et al 2005). Sensitive surveys at the latter wavelengths have only 
recently become possible thanks to the Spitzer satellite. 

For these reasons. Glass et al. (120011 1 conducted a /T-band 
(2.2 /im) survey for variable stars covering 24 x 24 arcmin^ (56 
X 56 pc at a distance of 8 kpc) and centered on the GC in a 
study spanning 4 years. The majority of the variable sources they 
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found were, as expected, Miras and OH/IR stars with periods 
ranging from 150 d to about 800 d. Uncertainty as to the fore- 
ground extinction unfortunately precluded any detailed compar- 
ison of their luminosities with similar populations in other well- 
studied areas, such as the the solar neighbourhood, the Baade's 
Windows and the Magellanic Clouds, where period-luminosity 
relations have been determined. 

The inner Galaxy has also been searched intensively in the 
radio region for OH sources (see e.g. Lindqvist et al. 119921 
Sjouwerman et al. ,1996. Wood et al. il998 Blommaert et al. 
II998I Vanhollebeke et al. 12006b . From these surveys, it is clear 
that a number of the large amplitude variables remained unde- 
tected in the Near-infrared, owing to the extremely high extinc- 
tion in some regions. Nevertheless, from OH and SiO observa- 
tions, (see e.g. Messineo 120041 and Deguchi et al. ,2008 ) radial 
velocities are available for many of these sources, making them 
extremely valuable for studies of stellar kinematics near the GC. 

Because the large-amplitude variables have very high mass- 
loss rates they also play an important role in the chemical evo- 
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lution of the Galaxy. It is therefore desirable that they should be 
characterised as fully as possible. 

Previous studies in the near-IR using 2MASS or DENIS data 
have yielded extinction maps with a typical resolution of several 
arcmin (see e.g. Schultheis et al. il999i Marshall et al. I2006I I. 
These maps, however, were limited by the sensitivities of the 
near-IR surveys and are only realistic in regions where Ay is 
less than 25 magnitudes. Unfortunately, this limit is generally 
exceeded towards the GC. 

In this paper we show that an improved map of the inter- 
stellar extinction around the GC can be made by fitting mid- 
IR data on AGB/RGB stars obtained by the IRAC camera on 
Spitzer to isochrones (see also Ganesh et al., in preparation). 
Extinctions have been determined for much more heavily ob- 
scured areas than was possible previously, especially towards 
the GC itself. We apply the new extinction values to a discus- 
sion of the long-period variables in the region, concentrating on 
the period-magnitude relations that they obey in the IRAC bands 
and their mass-loss rates. 

2. The data set 

2.1. The Spitzer IRAC Point Source Catalog of the Galactic 
Center (GALCEN) 

The central 2.0° x 1.4° of the GC have been mapped with 
Spitzer/IRAC between 3.6 and 8.0 /vm (Stolovy et al. 2006 and 
Stolovy et al. 2008, in preparation). Ramirez et al. (2008) per- 
formed point-source extraction on the IRAC data and published 
a confusion-limited catalogue of point sources that also included 
photometry from 2MASS. The IRAC magnitudes are refeiTed 
to as [3.6], [4.5], [5.8] and [8.0], corresponding to their central 
wavelengths. The average confusion limits are 12.4, 12.1, 11.7 
and 11.2 magnitudes, respectively, but can vary by 2 or 3 mag- 
nitudes within the survey. The whole catalogue (referred to here 
as GALCEN, see Ramirez et al. 2008 for more details) consists 
in total of about one million sources. Most of these show the 
characteristics of red giants or AGB stars, but there are several 
hundreds of extremely red sources which may be massive Young 
Stellar Objects (YSOs). A major problem that must be faced 
in the study of the LPV component is the bright-star limit of 
the IRAC camera. According to the Spitzer Observer's Manual 
(Version 7.1), saturation starts at 7.9, 7.4, 4.8 and 4.8 magni- 
tudes respectively. Saturated sources are however retained in the 
GALCEN catalogue but if a flux is greater than the limit it is 
flagged with the number '3' (see Ramirez et al. 2008 for a more 
detailed discussion). We have excluded photometry with flux 
flag values of 3 from our analysis. 

2.2. GALCEN vs GLIMPSE-II 

Recently, the GALCEN IRAC observations have 
been processed independently by the GLIMPSE-II 
team (see http://www.astro.wisc.edu/sirtf/docs.html). 
Different data reduction and source extraction proce- 
dures from those of the GALCEN survey were used (see 
http://www.astro.wisc.edu/glimpse and Benjamin et al. 2003). 
These data were kindly made available to us in advance of pub- 
lication by the GLIMPSE-II team. As discussed in the online 
document, they found systematic offsets between the GALCEN 
and GLIMPSE-II photometry that was larger than the combined 
uncertainty in both observations, at the bright and faint ends of 
the observed range. To look for any possible systematic effects 
on our work we have also calculated the extinction maps that 



we discuss in Sect. 3 using the GLIMPSE-II results. We find 
that the differences in the values of the extinctions derived 
from GALCEN and from GLIMPSE-II are smaller than the 
uncertainties in the determinations. 

The locations of the LPVs in the IRAC colour-colour dia- 
grams, the IRAC colour-magnitude diagrams, and the period- 
luminosity relations (see Sections 4—5) also differ slightly 
between the GALCEN and GLIMPSE-II photometry. The 
GALCEN and GLIMPSE-II photometry differ by -0.05 mag 
in the mean, with a standard deviation of 0.12 mag, over the 
brightness range of the LPVs. We find no statistically signifi- 
cant dependence of the difference between the GALCEN and 
GLIMPSE-II photometry on magnitude, or on the periods of the 
LPVs (which can be regarded as a proxy for their de-reddened 
magnitudes). The scatter in this difference, however, is large and 
we conservatively place a limit of 0.2 mag on any possible varia- 
tion in its value over the brightness range occupied by the LPVs. 
Our conclusions in this paper are, however, unaffected by photo- 
metric uncertainties at this level. 



2.3. Long-Period Variables 

The Glass et al. ( 1200 U survey covered an area of 24 x 24 
arcmin^ around the GC, as mentioned, and found about 400 pe- 
riodic variables with Mira-like amplitudes and an average pe- 
riod of 427 d. (Among these are 64 OH/IR stars included in the 
Wood et al. 1998 sample). We refer to the survey list hereafter as 
the Glass-LPV catalogue. Compared to the well-explored SgrI 
Baade's window of low extinction which lies at / - -1.4° and 
b = -2.6°, the GC field contains more than ten times the number 
of variables per arcmin-. The average period in the SgrI field is 
much lower, at 333 d. Due to the limited photometric precision 
of the Glass et al. (2001) survey, small amplitude variables such 
as semi-regulars were not detected. 

2.4. ISOGAL observations 

The details of the ISOGAL observational procedure with 
ISOCAM (Cesarsky et al. 1 19961 1 and the general processing of 
the data are described in Schuller et al. ( 12003b . In this paper we 
will discuss only the ISOGAL FC-00027 - 00006 field (here- 
after called FC-027), centered at (l,b) = (-0.27°, -0.06°). It 
was observed with the narrow filters LW5 (6.5-7. 1 fim) and LW9 
(13.9-15.9 fim), using 3" X 3" pixels, in order to avoid sat- 
uration. More specific details about this field can be found in 
Schuller et al. (.200 6 ). The parts of the region closest to the GC 
itself could not be observed by ISOGAL due to saturation prob- 
lems even when the naiTow-band filters were used. Nevertheless, 
25% of the Glass-LPV area is overlapped by the ISOGAL survey 
at 7 and 15 fim, mainly towards the south-western corner. 

2.5. Cross-identifications 

Cross-identification of stellar catalogues in these extremely 
crowded regions, which are often confusion limited, is a very 
difficult issue. In order to minimize spurious associations we 
used the following method when comparing the different cata- 
logues mentioned above. We first calculated the offsets of the 
400 brightest sources between a given pair of catalogues and 
determined the second-order distortion matrix. Following appli- 
cation of the offset and distortion corrections, we made the sub- 
sequent cross-identifications using a search radius of 2 arcsec: 
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Fig. 1. Upper panel: Interstellar extinction map of the central 280 x 200 pc of the Galactic Center with spatial resolution 2', derived 
using 2MASS J - data. Lower panel: As above, except using the [3.6]-[5.8] colour-magnitude data for point sources from the 
Ramirez et al. (2008) GC catalogue. 



First, we cross-identified the LPVs of Glass et al. (2001) 
with the GALCEN catalogue. By limiting the GALCEN cat- 
alogue to the area covered by the Glass-LPV survey, it was 
reduced to a total of 46314 sources. Of the 421 LPVs, 410 
were found in the GALCEN list (97%). Many of them are, 
however, saturated, especially at [3.6] and [4.5]. (From the 



Glass-LPV sample there were 162, 178, 36, and 58 saturated 
sources at [3.6], [4.5], [5.8], and [8.0], respectively). 
Second, we cross-identified the Glass-LPV catalogue with 
ISOGAL. We found 183 matches, despite the small overlap. 
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3. Interstellar extinction 

3.1. Derivation of the extinction vaiues 

As stated, interstellar extinction remains a serious obstacle to the 
interpretation of stellar populations in the GC region. 

The present work follows the method of Ganesh et al. (in 
preparation) who recently mapped the extinction towards the 
inner Bulge with a spatial resolution of 2'. GLIMPSE-II and 
2MASS data are compared to the latest isochrones for evolved 
stars (Marigo et al. I2008I I in order to derive Ay. Several studies 
of this type, using only the 2MASS data, have been undertaken 
before (e.g. Schultheis et al. [19991 Marshall et al. l2006b . 

In the method of Ganesh et al, colour-magnitude diagrams 
are constructed within sampling boxes of 2' square and the 
amount by which each individual data point has to be de- 
reddened for it to fall on the isochrone is determined. The ex- 
tinction values for Ay > 20 are derived using the [3.6]-[5.8] 
colour, while 2MASS J and data are used for Ay < 20. 

With the large number of filter combinations available (three 
from 2MASS and four from GLIMPSE-II) several possible com- 
binations had to be considered. The [5.8] vs. [3.6]-[5.8] colour- 
magnitude diagram was selected since it was found to yield the 
smallest dispersion, i.e. the [3.6] - [5.8] colour excess was found 
to be the most suitable one for determining the interstellar ex- 
tinction (see e.g. Indebetouw et al. 2005). This pair of filters has, 
in fact, the advantage of simultaneous observations (the other 
two IRAC colours were observed at different times) and also has 
better completeness (more stars detected through high extinc- 
tion). Further, at 3.6 yum and 5.8 //m, the observed flux is dom- 
inated by stellar photospheric emission while at 8.0;um inter- 
stellar PAH emission becomes prominent. The models indicate 
that metallicity effects are negligible in the [3.6]-[5.8] colour, in 
contrast to 7 - Ks, where the RGB/AGB branch of a metal-rich 
stellar population shifts to redder J - Ks colours (Schultheis et 
al. |2004). For more detailed information see Ganesh et al. (in 
preparation). 

In the following, we took the GALCEN catalogue of 
Ramirez et al. ( I2008I I and limited ourselves to stars with [3.6] 
< 12.5 and [5.8] < 12. We have also excluded sources with 
[3.6]-[4.5] > 0.8 and [5.8]-[8.0] > 0.5 from our extinction cal- 
culation, because sources with these colors are typical of YSOs 
(see Allen et al. 2004) and the inclusion of their extreme col- 
ors would bias the results. We again chose a sampling box of 
2' square in order to get a sufficient numbers of stars for the 
isochrone fitting. For each sampling box the average extinction 
value was taken. 

The values ofA^/Ay that we use to make the extinction maps 
were derived by using the infrared colour-colour diagram J - K, 
vs. /^^-[IRAC] as proposed by Jiang et al. ( 120031) . Indebetouw 
et al. ([2D05T l and Ganesh et al. (2008). We used the J - Ks vs. 
A'j-[3.6] and the J-Ks vs. /rj-[5.8] diagrams to determine A[3 g] 
and A[g 8]. We assumed that most of the sources are luminous 
RGB stars or AGB stars with moderate mass-loss and similar in- 
trinsic colours. We further adopted the values Aj = 0.256 X Ay 
and Aks - 0.089 x A ^ as used by Schultheis et al. (1999). Fitting 
straight lines in these diagrams gives us the following values: 
A[3.6]/Av = 0.0498 ± 0.0015 and A|5.8]/Av = 0.0308 ± 0.0015. 
These values were determined using the whole GALCEN cat- 
alog. Ganesh et al. (2008) will discuss more in detail how the 
extinction coefficients vary with different line of sights using the 
whole GLIMPSE-II dat set. Our A[3 5]/Av value agrees within 
the errors with those from Indebetouw et al. (2005), Lutz (1999), 
Flaherty et al. (I2007I I and Roman-Zufiiga et al. 120071 while our 



deiived A[5 gj / Ay is a litte bit lower, but still Ues inside the eiTors 
of the quoted authors. 

Lutz et al. (1996) and Lutz (1999) measured a flat mid-IR 
extinction curve towards the GC. However Flaherty et al. (2007)1 
and Roman-Ziiniga et al. 12007 l have shown that extinction coeffi- 
cents may vary somewhat depending in which line of sight one is 
looking (e.g. especially towards star forming regions). Recently, 
Chapman et al. (2008) have also studied the mid-infrared extinc- 
tion law in three molecular clouds and found that as the over- 
all extinction increases, the curve flattens out. The coefficients 
evidently vary according to whether the dust is located in the 
general ISM or in dense molecular clouds. 

Figure [T shows the Ay map derived from the 2MASS filters 
(upper panel) and that from the IRAC [3.6] and [5.8] filters of 
GALCEN (lower panel). There is a great similarity in the overall 
distribution of the two extinction maps. Where 2MASS gives 
lower limits for the extinction or is unable to probe the extinction 
due to lack of sources, the GALCEN extinction map exhibits 
high values for the extinction {Ay > 40). Within Figure T] there 
are regions of low extinction according to the 2MASS data but 
which are seen to have very high extinction when looked at by 
IRAC (e.g. the region around Z = 0.2 ° and b = -0.5 °). These 
areas are mostly associated with infrared dark clouds. 2MASS 
gives artificially low values of extinction in such cases because it 
detects only those stars located in front of the dark cloud whereas 
IRAC penetrates to much greater depth. 
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Fig. 2. Extinction values derived from 2MASS compared to 
those derived from IRAC colours. The straight line gives the 
identity relation. The dotted line indicates the maximum value 
of Ay to which 2MASS is sensitive. 



Since the [3.6] - [5.8] colour is less than J - K^, for a given 
amount of reddening the precision it yields in Ay is not as high. 
The loss of accuracy at low Ay can, however, be avoided by us- 
ing the near-infrared derived values for mapping when they are 
available. At larger Ay, the results from the [3.6] and [5.8] map 
are, of course, to be preferred. 
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Fig. 3. [5.8] vs. [3.6]-[5.8] colour-magnitude diagrams towards regions of low (left panel), intermediate (middle panel) and high 
(right panel) extinction (see text) in the GALCEN survey. The corresponding isochrones are displayed. The horizontal dotted line 
indicates the approximate confusion limit. The corresponding fields (from left to right) are a - 265.93596 °, 6 = -28.44771 °; a = 
266.16275 °, 6 ^ -28.44771 °; and a = 267.027222 °, 6 ^ -28.967888 °. The size of each field is 2'. The region with Ay = 70"' is 
located at a =267.027222 ° and (5 = -28.967888 ° and is associated with the infrared dark cloud IRDN3 of Dutra & Bica (2001) who 
list it as an opaque region at Kg with an angular diameter of 3'. 



Different ages or metallicities for the Bulge population do 
not significantly affect the mid-infrared colours of RGB stars, 
as will be seen. The sequences in the sampling fields are well 
defined in the case of low extinction while there is some scat- 
ter when it is high. This scatter arises from variable extinction 
on small scales within the Bulge and along the line of sight. 
Towards regions with very high extinction (Ay > 60), the num- 
ber densities in the CMD can drop significantly due to the loss 
of the fainter stars. Such points have relatively large errors. 

Figure |2] compares the extinction values derived from J and 
K, 2MASS data and those derived from the [3.6] - [5.8] IRAC 
colours. There is a tendency towards a linear relation between 
Av(2MASS) and Av(IRAC) up to Ay - 20 with a dispersion of 
Ay - 3. For higher Ay values, 2MASS can no longer detect the 
most reddened stars at the GC and is biased towards less red- 
dened foreground stars. For those regions where 2MASS gives 
low Ay values and IRAC high ones, different optical depths are 
obtained. The values of Ay derived from J - K, are representa- 
tive of the foreground population rather than the regions hidden 
by high extinction, while the IRAC bands relate to stars that are 
deeper within the clouds. 

Figure [3] shows the CMDs towards three regions of different 
extinction in the GALCEN survey. While for the low extinction 
field (left panel) the red giant branch is well defined and the dis- 
persion is low, the CMD is more scattered in the high extinction 
regions. Despite the high extinction, however, we get a sufficient 
number of sources to determine Ay. 

3.2. The combined extinction map 

Figure |4] shows the combined extinction (2MASS and IRAC) 
map as well as an intrinsic uncertainty map derived from the 
isochrone fitting. As explained before, for Ay{IRAC) < 20 the 



2MASS J - K, colour was used, while for Ay(IRAC) > 20 the 
IRAC [3.6]-[5.8] colour was used. The highest Ay values we 
can detect correspond to about Ay ^ 90. These are mostly asso- 
ciated with clouds that are dark even in the infrared. 

We will provide in electronic form (at CDS) the three extinc- 
tion maps as well as the error map in the form of individual FITS 
files. 

3.3. The uncertainty map 

The intrinsic uncertainty map corresponding to the combined 
interstellar extinction map is shown in Fig. |4] (lower panel). 
The uncertainty map is calculated as the standard deviation of 
the Ay distribution. For Ay < 20 the intrinsic uncertainty is 
about 2-3 mag while for Ay > 20 it can go up to 6-7 mag. 
Figure |5] shows the comparison between our Ay determinations 
based on GALCEN and on the GLIMPSE-II catalogue (see also 
Sect. 2.2). Despite the differences found in the photometry be- 
tween these two pipelines, no systematic offset is found between 
Ay derived from GALCEN and Ay derived from GLIMPSE- 
II. The typical r.m.s standard deviation between GALCEN and 
GLIMPSE-II is 6 mag in Ay which is the overall uncertainty in 
the derived extinction map. For very high Ay (Ay > 40), the 
differences between GALCEN and GLIMPSE become larger. 
Howver, for those regions the CMD is less well defined and the 
intrinsic uncertainty increases (see Fig.|4]). 

It should be noted that Nishiyama et al. (2006) found that 
the behaviour of near-IR extinction varies slightly though signif- 
icantly from one direction to another near the GC so that there is 
no universally valid extinction law for the region. These results 
were confirmed by Bandyopadhyay et al. (2008 ). Variations with 
direction of A^/Ay in the IRAC bands do not appear to be sig- 
nificant at the accuracy that is currently available, at least in the 
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Fig. 4. Upper panel: Combined interstellar extinction map of the central 280 x 200 pc of the Galactic Center. For AyilRAC) < 20 
the 2MASS J - colour was used, while for AyilRAC) > 20 the IRAC [3.6]-[5.8] colour was used. Lower panel: Intrinsic eiTor 
map of the combined interstellar extinction map. The eiTor is the sigma of the Ay distribution. 



fields studied by Indebetouw et al (2005) and Ganesh et al. (in 
prepration). 

Gosling et al. (2006) found complex small-scale structure in 
the infrared extinction towards the GC, with a typical size of 5- 
15", corresponding to 0.2-0.6 pc. They demonstrated that sig- 
nificant granularity appears to be present in high and variable 
extinction fields, indicating that extinction at J and H may have 
significantly higher values on smaller scales than those previ- 
ously derived for the GC as a whole. 



4. The Long Period Variables near the Galactic 
Center 

In the remainder of this work we discuss the long-period vari- 
ables of the Galactic Center region. Figure |7] shows the 2MASS 
image of the GC region at Ks- The known LPVs of Glass et 
al. (2001) are superimposed and are bright in Ks- We note that 
patchy and filamentary dark clouds (appearing white in Fig. 7) 
often vary on size scales smaller than the 2' sampling boxes used 
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Fig. 5. Extinction values derived from GALCEN compared to 
those derived from GLIMPSE-II. The straight line gives the 
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here to calculate the extinction. This can contribute to the pho- 
tometric errors in the extinction correction towards individual 
sources. 

4.1. Stellar isochrones 

The newest set of isochrones for the RGB/AGB phase include 
an improved treatment of the thermally pulsing asymptotic giant 
branch (TP-AGB) phase (Marigo et al. 2008 ). These isochrones 
predict realistic tracks for the TP-AGB, taking into account 
mass-loss and pulsation. The whole of TP-AGB evolution is now 
treated in a realistic way, especially the crucial effects of the third 
dredge-up, hot bottom burning and variable molecular opacities. 

The dust models incorporated in these isochrones are an ex- 
tension of the work by Groenewegen (12006b . They have been 
calculated with a 1 -dimensional dust radiative transfer code that 
solves the radiative transfer equation and the thermal balance 
equation in a self-consistent way. 

The main physical inputs to the model are the luminosity, 
distance, photospheric spectrum, mass-loss rate, dust-to-gas ra- 
tio, expansion velocity, dust condensation temperature and com- 
position of the dust. For oxygen-rich stars two main species of 
dust are considered: aluminium oxide (AI2O3) and silicate dust. 
Blommaert et al. (2006 ) studied the CVF ISOCAM spectra of a 
sample of AGB stars with low mass-loss in the Galactic bulge 
and found that the dust content is dominated by AI2O3 grains. 

Isochrones for different metallicities, Z=0.008, Z=0.019 and 
Z=0.038, and ages between 0.5 and 10 Gyr, have been calculated 
for the Spitzer IRAC filters. In the present work, we consider the 
isochrones placed at the distance to the GC (8.0 kpc, Reid ll993T l. 




([3.6]-[8.0])„ 

Fig. 6. ([3.6] - [8.0])o vs ([3.6] - [5.8])o diagram. Open circles 
indicate the GC LPVs of Glass et al. (2001). Filled red circles 
denote oxygen-rich LPVs in the LMC (see text). The straight 
lines are models computed by Groenewegen (2006) using realis- 
tic stellar atmosphere models. The dashed green line is a model 
with 100% aluminium oxide dust, the dotted blue line is 60% 
aluminium oxide and 40% silicate, and the solid magenta line 
is 100% silicate dust. The extinction vector of Ay - 20 is indi- 
cated. 



4.2. LPVs in the Spitzer colour-magnitude and colour-colour 
diagrams 

Figure [8] (left panel) shows the [3.6]-[8.0] vs [8.0] diagram. 
Indicated are also the approximate saturation limits at [3.6] and 
[8.0]. In total, only 248 known LPVs have reliable (unsaturated) 
IRAC photometry. On the right panel of Fig. [8] the dereddened 
[3.6]-[8.0] vs [8.0] diagram, derived using the extinctions given 
in Fig.|4] shows clearly that the locus of the known AGB stars is 
above the RGB tip ([8]o - 7.3). As can be seen, the isochrones of 
Marigo et al. (2008) are a good fit to the dereddened CMD. This 
agreement confirms the correctness of the derived extinctions. 
The locations of the isochrones change only marginally with age 
and/or metallicity (e.g. from 1 to 10 Gyr and 0.008 < Z < 0.03 
in these bands. The AGB stars show a large scatter in this dia- 
gram due to strong variability (we only have single-epoch mea- 
surements) and also to the infrared excesses arising from heavy 
mass-loss (see Sect. |6]l. A similar diagram can be drawn for 
LPVs at [4.5] where we find only 84 AGB variables to be un- 
saturated. 

4.2.1 . Comparison with Spitzer observations of LMC LPVs 

Cioni et al. (|2001) examined a field of about 0.5 square de- 
grees around the optical center of the Large Magellanic Cloud. 
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Fig. 7. 2MASS image (reversed grayscale) of the Galactic Center region at K^. The known LPVs from Glass et al. (2001) are 
superimposed with red circles. 
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Fig. 8. Left panel: [3.6]-[8.0] vs [8.0] colour-magnitude diagram of the GALCEN field, with the known LPVs (red crosses) and 
OH/IR stars (blue open circles) superimposed. The extinction vector for Av=20 is indicated. The dotted horizontal line indicates the 
saturation limit at [8.0] while the dotted diagonal line shows the saturation limit at [3.6]. Right panel: Same diagram but dereddened 
(see text). Superimposed is the isochrone for 10 Gyr and Z=0.02 put at a distance of 8 kpc. 



Light curves in the EROS V and R bands were obtained for 334 
variables, most of which are of long period (including Miras 
and semiregulars). In addition they obtained spectroscopy of 
about 100 AGB variables whose visible spectra were classi- 
fied as carbon- or oxygen-rich. We used their sample of known 
AGB stars and cross-identified their catalogue with the SAGE 



(Surveying the Agents of a Galaxy's Evolution) MIPS and IRAC 
imaging survey (Meixner et al. 2006). In the latter work, a field 
of 7 sq degrees of the LMC was surveyed between 3.6 and 160 
fim. In total, 33 1 objects were found in the SAGE catalogue. For 
further discussion, we will refer to this as the LMC-AGB sam- 
ple. 
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Figure |6] shows the dereddened [3.6]-[8.0] vs [3.6]-[5.8] 
colour-colour diagram for both the Galactic Center and the LMC 
AGB LPVs. We have superimposed for comparison theoretical 
dust models by Groenewegen ( 12006b . These utilise several types 
of dust for the O-rich AGB stars: 100% aluminium oxide, a com- 
bination of 60% aluminium oxide and 40% silicate, and pure 
silicate dust (see Groenewegen 2006 for a discussion). All three 
models fail to match the observed colours and a significant offset 
is clearly visible. 

Our comparison is restricted to oxygen-rich AGB stars (283 
stars). It should be noted that we split the LMC-AGB sample 
into oxygen-rich and carbon-rich AGB stars using only the J-K 
colour. Though Cioni et al. (.2001 ) showed that carbon stars can 
in general be eliminated by avoiding colours in the range 1 .4 < 
J - Ks < 2.0, they warn (see their Fig. 3) that stars selected in 
this way will contain a certain fraction of carbon stars. 

Figure |6] shows that the GALCEN variables and LMC 
oxygen-rich variables only partly overlap in the colour-colour 
diagram. The GALCEN LPVs show much redder colours in 
[3.6]-[8.0] and [3.6]-[5.8] than the LMC sample. These stars 
also have a higher proportion of long periods (Groenewegen & 
Blommaert 2005) relative to the LMC and the remainder of the 
Galactic Bulge. The redder colours in the IRAC bands may in- 
dicate that mass-loss rates of AGB stars in the GC are higher in 
general than in the LMC. However, the differing circumstellar 
dust-to-gas ratios between the LMC and GC makes direct com- 
parison difficult. The differences can also be a manifestation of 
the strong relation between period and mass-loss (see Sect.|6]l. 



5. The Period-magnitude relations 

AGB variable stars in the Galactic Bulge follow a well-defined 
period-magnitude relation in the near infrared. Glass et al. 
(2001) showed that the mean K magnitude for the Glass-LPVs, 
measured over a complete cycle, showed some tendency to fol- 
low a log P\s. K relation in spite of a very high dispersion due to 
the large and unknown variations in Ay. Some of the scatter may 
also have arisen because of crowding effects or insufficient ob- 
servations having led to incorrect periods, especially at the long 
(log P > 2.8) and short ends (log P < 2.3) of the range. The 
Glass et al. (2001) sample contains a higher fraction of long pe- 
riods than e.g. the SgrI field. Thus, circumstellar reddening due 
to heavy mass-loss affects the log P vs. K relation by depressing 
the emergent fluxes at longer wavelengths. 

Figure [9] shows the IRAC magnitude vs. log P relations 
for the Glass-LPVs, after correction of the photometry for 
interstellar extinction. We find following relations for the 
Glass-LPV sample: 
[3.6] = log Px -5.21 -1-18.94 
[4.5] = log Px -5.71 4-19.98 
[5.8] = log Px -6.10 -1-20.70 
[8.0] = log P X -6.73 -H21.96 

The typical r.m.s deviation for the four [IRAC] bands are 
0.71 , 0.62, 0.60 and 0.61 mag in the period range between 
2.2 < logP < 2.8. This is smaller than at Kq where the rm.s 
is 0.95 mag (corrected according to the extinctions determined 
in this work) as the IRAC bands are less affected by any residual 
errors in the extinction. Even within the IRAC bands, the disper- 
sions at the longer wavelengths are less than that at 3.6 fim. The 
Glass-LPVs with periods greater than 900 days do not follow the 
log P vs. [IRAC] relations. Their periods could be incorrect as 
these only have a few observations (see Glass et al. 2001) and in 



addition these stars suffer heavy mass-loss which depress their 
[IRAC] fluxes. 

Superimposed on Fig. |9] is the LMC-AGB star sample. 
Again, we have chosen only oxygen-rich Mira variables in or- 
der to make direct comparisons with the Glass-LPV sample. We 
have put the LMC LPVs at the distance of the GC using a dis- 
tance modulus difference of Am = 4.0 mag. Within the errors, 
the Glass-LPVs do follow the LMC relation. 

As mentioned, the scatter in the period-magnitude relations 
for the Glass-LPVs is much higher than in the LMC-AGB sam- 
ple. This is due to the uncertainties in the extinction determina- 
tion. No attempt has been made to compensate for the triaxial 
distribution of stars in the Bulge. Recent studies of the central 
region of our Galaxy suggest a smaller bar-like structure within 
the main bar (see e.g. Alard 2001 , Nishiyama et al.'2005) which 
might be partly responsible for the observed scatter However, 
the parameters of the central bar (such as its angle, length and 
orientation) are still poorly known. 

In contrast to the light curves obtained by Glass et al. (2001), 
we have only single epoch measurements in the IRAC bands. 
The K amplitude distribution of the Glass-LPVs peaks around 
0.8 mag, going up to 1.5 mag for the OH/IR stars. In order to es- 
timate the expected amplitude of LPVs in the IRAC filter bands, 
we used the calculated colours of Galactic oxygen-rich AGB 
stars by Marengo et al. (2006). Several such objects were ob- 
served by ISO at different epochs. Mira variables, such as Z Cyg, 
can show variations up to 1 mag in each IRAC channel, while in 
the most extreme cases, such as the dust-enshrouded OH/IR star 
VI 300 Aql, variations of 1.4 mag may occur in each filter 

Up to now, period-magnitude relations have been mainly 
studied in the intermediate Bulge (such as Baade's Window) 
and the LMC, where the log P vs Ks relations seem to be in- 
dependent of the galactic environment (see e.g. Schultheis et al. 
2004, Glass & van Leeuwen 2007). Groenewegen & Blommaert 
(2005), by using OGLE data, also compared the period- 
magnitude relations in the Galactic bulge and the Magellanic 
Clouds and concluded that there is no difference in their slopes. 
Rejkuba 620041) found that even Miras in Cen A follow a similar 
relation. 

The Galactic Center, however, has up to now been excluded 
from this comparison. In asking whether substantially different 
behaviour should be expected we first summarize recent abun- 
dance studies in order to place the GC in context with the Bulge 
([Fe/H] ^ -0.1; alpha-enhanced) and the LMC ([Fe/H] ^ -0.4; 
no alpha-enhancement). Studies have shown that in the GC the 
mean [Fe/H] is nearly solar; [Fe/H] ^0.1 (Ramirez et al. 120001 
Najarro et al. 2004 2008 ; Cunha et al. 2007 ). On the other hand, 
Ramirez et al. (2000) found that the distribution in metallicity is 
significantly narrower than in the Bulge, emphasizing the differ- 
ence between its stellar population and that of the GC. 

Figure |9] shows clearly that the slopes of the period- 
magnitude relations at the IRAC wavelengths for the Glass- 
LPVs are similar to those of the LMC. This suggests that any 
dependence of the log P vs. IRAC relations on abundance, if 
present at all, must be small. Within the uncertainties, there 
is no evidence for a significant difference between the period- 
magnitude relations in the LMC and in the GC. This is in agree- 
ment with Whitelock et al. ( 12008t who found, after reanalyz- 
ing published lightcurves of AGB variables in the LMC, similar 
zero points in the period-Mj^ relations for systems with differ- 
ent metallicities. They did not, however, include in their analysis 
the variables towards the GC, which is thought to be the most 
metal-rich environment. Thus, our analysis shows that the appar- 
ent universality of the PL relations extends to the IRAC bands. 
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Fig. 9. Dereddened IRAC magnitudes (using the extinction map as described in Sect. 3) vs. log P for LPVs in the Glass et al. (2001) 
GC field (Glass-LPVs; black dots), where P is the period. Open magenta squares are LMC-AGB stars (oxygen-rich Mira variables). 
The straight red line is a least-squares fit to the LMC-AGB data set. 



6. Mass-loss 

AGB stars contribute more than 70% towards the enrichment of 
the dust component of the interstellar medium (ISM) in the so- 
lar neighbourhood (Sedlmavr Il9941 l and hence it is important 
to study their mass-loss in other parts of the Galaxy. One of 
the most promising tools for determining mass-loss rates is the 
combination of near-IR and mid-IR colours such as the IRAS 



Kq - [12] or the ISOGAL Kq - [15] colours (see e.g. Whitelock 
et al. 1994; Habing 1996; Le Bertre & Winters 19 98; Omont et 
al. 1999; Jeong et al. 2002, Ojha et aL |2003ll2007T l. 



We used the dust radiative transfer code for oxygen-rich 
AGB stars from Groenewegen (2006) to derive theoretical mass- 
loss rates. The procedures described by Ojha et al. (120071 ) to de- 
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rive observed mass-loss rates based on the (Ks 
(see Equ. 1). 



[15])o colour 



log M 



-8.6171 + 0.85562 X -0.064143 + 0.0018083x^(1) 



where X = (/:- [15])o 

As in Ojha et al. (12007!) we are considering an AGB star 
model with T^ff = 2500 K and 100% silicate composition. Note 
that the empirical values determined by Groenewegen (2006) are 
for mass-loss rates smaller than 2 x 10"^ M© yr"'. For the high 
mass-loss end we extrapolate the empirical relation. 

The distribution of mass-loss rates is displayed in Fig. [TO] 
The derived values of M range from 10"^ to IQ-^ Mq yr '. The 
distribution, however, is definitely incomplete for M < 1 x 10"^ 
Mq yr"'. We emphasize that the determination of these rates is 
very uncertain as discussed by Schultheis et al. (120031) . The un- 
certainties could easily reach a factor of 3. 
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Fig. 10. Distribution of mass-loss rates for the GC Glass- 
LPV sample. The mass-loss rates are inferred from models 
(Groenewegen 2006) using the dereddened A'- [15] colour 

There is some correlation between mass-loss rates and 
period, in the sense that longer periods show higher rates. 
Theoretical isochrones by Marigo et al. (2008) which include 
mass-loss and pulsation periods predict that there should be such 
a correlation; however there is still some disagreement with the 
observational data indicating that self-consistent dust models of 
oxygen-rich AGB stars are in need of further improvement. 

7. Luminosities 

We have derived a bolometric magnitude (Mboi) for each Glass- 
LPV from its de-reddened K magnitude, an assumed bolometric 
correction (BCa:), and a distance modulus of 14.5 for the GC. 




Fig. 11. Histogram of luminosities of the GC Glass-LPV sample 
(solid line). The luminosities for an inner Bulge sample (Ojha 
et al. 2007) is shown for comparison (dotted lines). The vertical 
lines mark the tip of the RGB (at fainter Mboi) and the tip of the 
AGB (at brighter Mboi). 



We used an approximate relation between BCk and (K-[l5])o 
derived by Ojha et al. (2003 ) for AGB stars, using data from 
Frogel & Whitford (1987). We refer here to Ojha et al. (120031) for 
a more detailed description of the method. These authors con- 
clude that the rms uncertainty in the luminosities of the Bulge 
sources that they observed with ISO is about 0.5 mag. 

Fig. [TT| shows the luminosity function of the Glass-LPVs. 
These include high luminosity stars, with Mboi < -6, which are 
not observed in the inner Bulge AGB stars (Ojha et al. 2007). 
However, one should note that the sample of Ojha et al. (2007) 
excludes the GC region and are only single-epoch measure- 
ments. The luminosity distribution of the Glass-LPVs, compared 
to that of the inner Bulge AGB stars (Ojha et al. 2007), is also 
broader and less strongly peaked. Due to the variability of our 
sources, the luminosity distribution will be widened. Thus, com- 
pared to the Bulge sample of Ojha et al., we find an excess in the 
GC of luminous stars. 

Such an excess of luminous stars in the GC, compared to 
Bulge fields such as Baade's Window, has previously been found 
in the infrared stellar luminosity function by Blum et al. (. 19961 ). 
Narayanan et al. (119961) . Launhardt et al. (12002b . Genzel et al. 
(2003ll, Van Loon et al. (^03 ), and Figer et al. (l2004l l. In ad- 
dition, Blommaert et al. (1998) and Wood et al. (1998) found a 
similar excess of high luminosity OH/IR stars in the GC com- 
pared to the Bulge. All these authors conclude that this is gen- 
erally evidence for recent star formation in the GC. Blum et al. 
(11996 ). Figer et al. (2004 ), and Maness et al. (2007), among oth- 
ers, have derived detailed star formation histories for the GC by 
fitting models to the observations. Their studies suggest that star 
formation has occurred within the last 10 Myr within the cen- 
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tral 30 pc, and that it has been more or less continuous over the 
Hubble time. 

Groenewegen & Blommaert ( 12005b studied the Mira period 
distribution of six fields at similar longitudes but spanning lat- 
itudes -5.8 < b < -1.2°. They found for all these fields a 
similar period distribution, consistent with the average period 
in Baade's Window of ~333 d (Glass et al. 1995 ). By comparing 
their fields with the LPVs of Glass et al. (2001 ) they found that 
there is a significant overabundance of LPVs with periods ^500 
d in the GC compared to Bulge fields at higher latitudes. This 
excess of LPVs with long periods in the GC is not surprising, 
considering the excess of luminous sources in the GC generally 
(see Fig. 
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when combined with the period-luminosity rela- 
tion. Groenewegen & Blommaert (12005 l l explain this difference 
in the period distribution as due to a younger population in the 
GC with a higher initial mass in the range of 2.5-3 Mq compared 
to 1-2 Mq for higher latitude Bulge fields. 

We conclude that the GALCEN-LPVs show an excess of lu- 
minous, young stars compared to the galactic Bulge. It is un- 
fortunate that the isochrones in the [IRAC] bands are not very 
sensitive to metallicity and/or age; thus no estimation of initial 
masses or metallicities in this wavelength range can be made. 



8. Conclusions 

We have presented a high resolution interstellar extinction map 
of the Galactic Center (GC), derived from the Spitzer IRAC 
catalogue of GC point sources (Ramirez et al. 2008). It com- 
bines observations of the RGB/AGB population with the newest 
isochrones (Marigo et al. 2008). We have improved on previ- 
ous near-lR extinction maps by using extinctions derived from 
mid-IR colours to resolve all regions of high Ay. The maximum 
value we found is Ay ^ 90. By using this extinction map, we 
have studied the Spitzer (IRAC) properties of the long-period 
variables (LPVs) of in the GC. They follow well-defined period- 
magnitude relations in the IRAC bands and are similar to those 
observed in the LMC, lending support to the suggestion that the 
log P vs. [IRAC] relations are universal, independent of metal- 
licity. Further, mass-loss rates for the Glass-LPVs have been de- 
termined by using additional data from ISOGAL. These LPVs 
show some correlation between mass-loss and pulsation pe- 
riod as predicted by theoretical isochrones (Marigo et al. 2008). 
Finally, the long periods of the Glass-LPVs as well as their lu- 
minosity function agree with the suggestion of Groenewegen & 
Blommaert (2005) that one is dealing with a young and more 
massive population in the GC. 
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